Crystal potentials under invariant periodic boundary conditions at infinity 
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The definiteness of bulk electrostatic potentials in solids under periodic boundary conditions defined 
in an invariant manner has been proved in the general case of triclinic symmetry. Some principal 
consequences following from the universal potential correction arising are discussed briefly. 



PACS numbers: 41.20.Cv, 61.50.Lt 
I. INTRODUCTION 

The summation of Coulomb potentials over crystal lat- 
tices is a classical problem Q important for determining 
the cohesive energy in crystals J2] as well as for describing 
electronic properties there Despite a long history 

of this subject (||] and references therein) the question 
about the uniqueness of electrostatic properties in the 
bulk remains controversial. Indeed, it is well-known that 
every electrostatic task is defied by boundary conditions 
and periodic boundary conditions are appropriate to 
the solution in the bulk 0. As a result, any electric field 
generated by polar unit cells is to be excluded as irrele- 
vant to the bulk state Notwithstanding, as far as 
electrostatic potentials are concerned, an arbitrary con- 
stant potential could be formally added as an available 
periodic solution of Laplace's equation 10 -1^]. There- 
fore, the absence of a constant potential in the original 
approach of Ewald is often treated as optional 
whereas the wide-spread standpoint is that the electro- 
static potentials are well defined by boundary conditions 
on open surfaces of crystals, but these potentials are un- 
determined to an additive constant in infinite crystals 
[Q. This claim is, however, at variance with at least 
the two physical statements: On the one hand, due to 
statistical arguments, the bulk state cannot be governed 
by surface ones |]l5| , ^ . Furthermore, the actual charge 
distribution in the bulk of a crystal is to be in one-to-one 
correspondence with potentials it generates 0,|8| . 

In the present paper we show that periodic boundary 
conditions defined properly arc sufficient to make elec- 
trostatic potentials definite in the bulk. It means that 
periodic boundary conditions are as effective as bound- 
ary conditions on surfaces are. Moreover, the latter ones 
can then be reconstructed as relevant to real surfaces 
consistent with bulk states [0Jl!I,E(]|l . 



II. 



LATTICE SUMMATION AS A CONVERGENT 
PROCEDURE 



In the general case of triclinic symmetry, let a crystal 
be described by primitive translation vectors a = ae^, 
b = hei, and c = cSc, where Ga, e;, and Ec are the ap- 
propriate unit vectors, with the products {ea^b) = cos a, 
(ebSc) = cos (3 and {gc^o) = cos 7. In terms of a charge 
distribution p(r) contained in a imit cell parallelepiped. 



the direct lattice-sum contribution to the electrostatic 
potential at a reference point r can be written as 



, Jv 



p(r') dr' 



|R,; 



- r — r 



(1) 



where i runs over the corresponding Bravais lattice spec- 
ified by Ri, the prime on the summation sign implies 
missing the singular contributions of the summand, the 
integration is over the unit-cell volume. 

To make the result of summation in (|^) definite, the 
following conditions for the absolute convergence of (|^) 
are to be suggested pl|,p2[: 



/ p(r) dr = 0, (2) 
Jv 

M^= [ r^p(r) dr = 0, (3) 
Jv 

Gf,i, = / r^r-yp(r) rfr = at n ^ ly, (4) 
Jv 



(5) 



where are Cartesian components of r, H is an arbitrary 
constant. If any initial p'"'(r) is not subjected to 
provided that the neutrality condition (j^) holds, then we 
may modify p"^'(r) as follows: 



p{r) 



(6) 



where S{r) is the Dirac delta function, qj are some ficti- 
tious point charges located in the vicinity of the 
origin, r = 0, at points rj connected by admissible lattice 
translations. The values of qj are constrained by 



(7) 



so as to exclude the contribution of qj to the overall initial 
charge distribution upon combining p(r) attributed to 
neighbouring unit cells. It is significant that ten different 
charge species among qj are sufficient to fulfil and 
(0). Thus, H remains optional in (||). 

It is advantageous to consider a unit cell with the origin 
in its geometric centre. Then a compact distribution of 
qj is supplied by Vj belonging to the following set of vec- 
tors: 0, ±a, ±b, ±c, ±a±b, ±b±c and ±c±a. Keeping 
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in mind that /3'"'(r) is still contained in the unit-cell par- 
allelepiped, p(r) can be connected with a parallelepiped 
that is the same in shape, but twice as large in size. It 
is convenient to adopt this parallelepiped as a new unit 
cell. Being additive, the potential of interest is then de- 
scribed by the sum of independent contributions of eight 
interpenetrating lattices composed of new unit cells each. 
Thus, it is conceptually sufficient to discuss the potential 
effect generated by a single lattice specified by (^-(1]), 
provided that this lattice is still determined by the lattice 
parameters a, b and c, as is verified later on. 



III. INVARIANCE OF PERIODIC BOUNDARY 
CONDITIONS 

To make the solution of interest determinate, we as- 
sume that the overall structure is composed of an integral 
number of unit cells restricted by planes which are paral- 
lel to the unit-cell faces and are specified by the vectors 
±A, ±B and ±C relative to a central unit cell, where 
A = Asa, B — Beb and C — Cbc at 



A B C ^ 

- = y = - » 1, 
a c 



(8) 



so that the uniformity along each crystallographic di- 
rection is maintained. Periodic boundary conditions are 
then readily involved as imposed in such a way that each 
couple of remote parallel restricting planes merges, so 
that equal number of complete unit cells occur along each 
direction of e^, e;, and Be. In this event, the invariant 
character of periodic boundary conditions implies that 
each plane of merging may also occur somewhere within 
boundary unit cells, without changing the result. It is 
important that if planes of merging happen in interme- 
diate positions specified by < / < 1 within boundary 
unit cells, without loss of generality, some instantaneous 
charge distributions are to be introduced on those planes 
so as to fulfil conditions (||) and (||) furnishing the conver- 
gence of the surface potential contributions. As a result, 
the contribution of any dipolar polarization along those 
planes can be eliminated, but dipolar moments normal 
to the planes in question are inevitable and contribute 
to the potential value in the interior 1 19 2^-|2^ . In the 
particular case of the -l-A plane one can show that the 
corresponding potential contribution takes the form 



<?a(/) = 



sin/3 



t(/) 

dt dp dup{t,p,u) 

-a/2 J-b/2 J-c/2 



ie{A} 



where in the triclinic co-ordinates r = tea + P^b + uBc, 
t{f) = af — a/2, i runs over unit cells truncated by the 
boundary plane, with R:t'^, the component of Ri along 
an outward normal to this plane. 



f? = [l — cos^ a — cos^ f3 — cos^ 7 



2 cos a cos P cos 7] 



1/2 



(10) 



Upon merging the ±A planes, auxiliary charges on those 
planes cancel each other at a given /, but the aforemen- 
tioned invariance of the boundary conditions is based on 
(^) averaged over /, so that the effective potential gen- 
erated by that couple of planes is to be specified as 



<Z'a(/)+<?-a(1-/) 



df, 



(11) 



where follows from (|^) upon inversion of the co- 

ordinate system. Carrying out the integration over / in 
(|n|) and utihzing (|), we obtain 



sin/3 



E 



R 



an Rl 



(12) 



IV. UNIQUENESS OF BULK POTENTIALS 

The structural factor described by the sum in (|l^) de- 
termines the potential at a large distance from the plane 
at hand. Therefore, it is independent of the discrete char- 
acter of that sum |ll|,|l9|,||j25) and so can be represented 
in the following integral form 



E 

^e{A} 



Ri 



if2 



Rf 



6c sin /3 7- 



dp 



du 



-c W^{a,p,uy 



(13) 



where the scale transformation to the parameters of a 
unit cell is performed in agreement with (0), 



W{a, 6, c) = (a^ 5^ + + 2ah cos a + 26c cos /3 



2ca cos 7 



1/2 



(14) 



Carrying out the integration in (|3|) and substituting the 
result into (O), we derive 



H 

V 



y(a,6,c|a,/3,7)-y(-a,6,c|a,/3,7) 
-r(a,-6,c|a,/3,7)-y(a,6,-c|a,/3,7) 
where v — abcf2 is the volume of the unit cell, 
y(a, 6, c|ci!,/3,7) = tan^^l [6csin^/3 a6r(7) 
-f^ car(a) - a^r{f3)] [a fl W(a, 6, c)] "^|, 

P(01 ) = cos </>! — cos (/)2 COS , 



(15) 



(16) 
(17) 



the parameters are the angles a, /3 and 7 in an arbi- 
trary combination. 

As a generalization of (p^, the total potential contri- 
bution associated with the ±A, ±B and ±C boundary 
planes after their merging takes the form 



2 



top 



(18) 



where <?b and <Pc are obtained from (|ll)-(0) upon the 
cychc interchanges A— »B^C, a^6—>c—>a 
and a — > /3 — > 7 ^ a there. To proceed further, we 
remark that the angles Y{a, b, c\a, (3, 7), Y{b, c, a|/3, 7, a) 
and y (c, a, 6I7, a, /?) associated with Vl^(+a, +&, +c) and 
denoted as Xj form a closed set with the property 

tan(xi +X2)tanx3 = 1, IXjl < 7i"/2. (19) 
One can readily prove therefrom that 

(20) 



= Xl + X2 + X3 



TT 

2' 



where the sign combination specifying the arguments 
of W{+a,+b,+c) is indicated as a superscript. Like- 
wise, for the angles associated with W^(— a, +6, -l-c), 
W^(+a, —6, +c) and W^(-fa, +6, -c) we get 



X 



X 



X 



TT 

2' 



(21) 



Substituting ( |20D and ( |2l| ) into (^, we finally reach 

27riJ 



<?top = 



(22) 



On combining (0) and (^2|), the definite bulk potential 
field takes the form 



C/b(r) -C/cd (r )+<?'. 



top ■ 



(23) 



Following Bethe |2^], one can see that the value of f7cd(r) 
averaged over a unit cell cancels the last term on the 
right-hand side of (p3|), so that for the mean bulk poten- 
tial we obtain 



C/b = 0. 



(24) 



Hence, the bulk potential field f7b(r) is independent of an 
optional parameter H and has no uniform component, in 
agreement with the result of Ewald M . 



V. DISCUSSION 

It is significant that in a general case of (^) the ef- 
fect of eight sublattices mentioned above on the result- 
ing <Ptop just compensates the increase of the unit-cell 
volume in each of them, so that relation ( p^ ) is repro- 
duced with the parameters attributed to the initial unit 
cell. It is also clear that the auxiliary fictitious charges 
introduced in (^) vanish under periodic boundary condi- 
tions, so that real structural charges are of importance 
altogether. On the other hand, according to the above 
procedure of averaging, the issue (^ turns out to be in- 
different to the particular definition of a unit cell. As dis- 
cussed in ^ , the foregoing result may also be associated 
with the translational invariance as an integral property 



as 



of the direct lattice sum (|^) within a special mode of sum 
mation. This circumstance was stressed by Ewald 
desirable upon the definition of lattice sums as such. 

Note that the potential C/cd(r) as a function of r is 
asymmetric if iJ 7^ 0. This is especially prominent in 
diatomic structures composed of point charges |^ , ^ , 
but relation (|3|) retrieves the symmetric result there. 

For completeness, one can show that, in terms of (|2 
the bulk Coulomb energy per unit cell takes the form 



(25) 



which is invariant, though p'"'(r) occupying a volume 
^ini j.g]jiains optional. According to [||, the variational 
derivative of 5b with respect to p'"'(r) is then equal to 



(5p""(r) 



t^b(r) 



(26) 



that is the basic statement for determining p'"'(r) in a 
self-consistent manner [ p^ . 

It is worth noting that relations ( p^ ) and (^3|) also 
describe the potential effect exerted by extended charges 
even if we deal with spherical electronic distributions in 
ions 1^. As a result, the potential asymmetry beyond 
the potential contribution of a point-charge lattice takes 
place as well. According to (|2^), this fact results in the 
asymmetry of the concentration of vacancies of different 
ionic species and so explains the n-type conductivity 
in intrinsic semiconductors such as ZnO or GaAs |pO|,pi| . 

Note that the description based on (||)-(0) and (|22|)- 
( p6| ) is quite general, with including the effect of the 
Lorentz field for polar unit cells as a particular case . 
A subtle problem associated with the definition of a lo- 
cal polarization in ferroelectrics |33|j3^ ] can also be eluci- 
dated therefrom, as will be discussed elsewhere. 

VI. CONCLUSION 

Without loss of generality, the problem of summation 
of Coulomb potentials over crystal lattices is investigated 
in terms of absolutely convergent sums with an arbitrary 
choice of the charge distribution in a unit cell. The prin- 
cipal case of triclinic symmetry is considered. It is shown 
that periodic boundary conditions imposed in an invari- 
ant manner so as to exclude the influence of the partic- 
ular choice of a unit cell are sufficient for determining 
the electrostatic potentials in the bulk as uniquely de- 
fined, with zero mean bulk potential value. A few direct 
consequences of the results obtained are pointed out. 

ACKNOWLEDGEMENTS 

I am grateful to Professor V. L. Ginzburg, Professor 
E. G. Maksimov and Professor V. G. Vaks for their en- 
couragement of this work. 



3 



* El-mail: kholopov@caspcr.che.nsk.su 
[1] M. L. Glasser and I. J. Zucker, in: Theoretical Chemistry: 
Advances and Perspectives, edited by H. Eyring and D. 
Henderson (Academic Press, New York, 1980), Vol. 5, 
pp. 67-139. 

[2] M. P. Tosi, in: Solid State Physics, edited by F. Scitz and 
D. TurnbuU (Academic Press, New York, 1964) Vol. 16, 
pp. 1-120. 

[3] F. E. Harris, in: Theoretical Chemistry: Advances and 
Perspectives, edited by H. Eyring and D. Henderson 
(Academic Press, New York, 1975), Vol. 1, pp. 147-218. 

[4] J. Ihm, Rep. Prog. Phys. 51 (1988) 105. 

[5] E. V. Kholopov, Preprint 2001-01 (Inst, of Inorg. Chem., 
Novosibirsk, 2001) (submitted to Philos. Mag. B). 

[6] J. D. Jackson, Classical Electrodynamics (John Wiley 
and Sons, New York, 1962). 

[7] M. Born and Th. von Karman, Phys. Z. 8 (1912) 297. 

[8] E. R. Smith, Proc. R. Soc. London A 375 (1981) 475. 

[9] M. W. Deem, J. M. Ncwsam and S. K. Sinha, J. Phys. 
Chem. 94 (1990) 8356. 
[10] J. P. Dahl, J. Phys. Chem. Solids 26 (1965) 33. 
[11] A. Redlack and J. Grindlay, J. Phys. Chem. Solids 36 
(1975) 73. 

[12] X.-G. Zhang, W. H. Butler, J. M. MacLaren and J. van 

Ek, Phys. Rev. B 49 (1994) 13383. 
[13] P. P. Ewald, Ann. Phys. F4 64 (1921) 253. 
[14] L. Kleinman, Phys. Rev. B 24 (1981) 7412. 
[15] M. Born and K. Huang, Dynamical Theory of Crystal 

Lattices (Clarendon Press, Oxford, 1954). 



[16] R. B. Griffiths, Phys. Rev. 176 (1968) 655. 

[17] M. A. Spackman and R. F. Stewart, in: Chemical Appli- 
cations of Atomic and Molecular Electrostatic Potentials, 
edited by P. Politzer and D. G. Truhlax (Plenum Press, 
New York, 1981), pp. 407-425. 

[18] K. A. Van Goncchton, W. J. Mortier and P. GeerUngs, J. 
Chem. Phys. 86 (1987) 5063. 

[19] P. W. Tasker, J. Phys. C 12 (1979) 4977. 

[20] D. Wolf, Phys. Rev. Lett. 68 (1992) 3315. 

[21] H. M. Evjcn, Phys. Rev. 39 (1932) 675. 

[22] C. K. Coogan, Aust. J. Chem. 20 (1967) 2551. 

[23] F. C. Frank, Philos. Mag. 41 (1950) 1287. 

[24] W. W. Lee and S.-L Choi, J. Chem. Phys. 72 (1980) 6164. 

[25] D. M. Heyes and F. van Swol, J. Chem. Phys. 75 (1981) 
5051. 

[26] H. Cokcr, J. Phys. Chem. 87 (1983) 2512. 

[27] H. Bethe, Ann. Phys. F4 87 (1928) 55. 

[28] L G. Fisenko and E. V. Kholopov, Phys. Status Solidi B 

173 (1992) 515. 
[29] E. V. Kholopov, Zh. Struc. Cliim. (to be published). 
[30] W. Van Gool and A. G. Piken, J. Mater. Sci. 4 (1969) 

105. 

[31] W. A. Gates, G. Eriksson and H. Wenzl, J. Alloys Comp. 
220 (1995) 48. 

[32] E. V. Kholopov, Preprint 2001-02 (Inst, of Inorg. Chem., 
Novosibirsk, 2001) (submitted to Philos. Mag. B). 

[33] D. Vanderbilt and R. D. King-Smith, Phys. Rev. B 48 
(1993) 4442. 

[34] G. Ortiz and R. M. Martin, Phys. Rev. B 49 (1994) 14202. 



4 



